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Abstract 

The cosmological baryon asymmetry can be explained by the nonperturbative electroweak 
reprocessing of a lepton asymmetry generated in the out-of-equilibrium decay of heavy 
right-handed Major ana neutrinos. We analyze this mechanism in detail in the framework 
of a SO(10)-subgroup. We take three right-handed neutrinos into account and discuss 
physical neutrino mass matrices. 
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1 Introduction 



One of the most striking features of the observable universe is the baryon asymmetry, which 
is usually expressed as the ratio of the baryon density % to the entropy density s. From 
measurements of the abundances of the light elements one findsQ: 

Y B = V*. = (0.6 - 1) x 1(T 10 . (1) 

s 

With the appearance of grand unified theories (GUTs) it became possible to explain this asym- 
metry by the baryon number (B) violating decays of Higgs or gauge bosons at the GUT scale. 
However, these models of baryogenesis are not easily reconciled with inflation. Indeed, a baryon 
asymmetry present before inflation would be diluted by a huge factor, while the reheating tem- 
perature after the inflationary phase is in general too low for these baryogenesis mechanisms 
to work. 

Then, it was realized that anomalous baryon number violating processes are unsuppressed at 
high temperatures 0j. These so called sphaleron transitions violate (B + L) and conserve (B — 
L), where L is the lepton number. As sphalerons are in thermal equilibrium for temperatures 
between ~10 12 GeV and ~10 2 GeV, they will strongly modify any primordial (B+L) asymmetry. 
The connection between the baryon asymmetry and a primordial (B — L) asymmetry is given 
by i 

Y - ( 8N f + 4NH ) Y (2) 

where Nf is the number of fermion families and Nh is the number of Higgs doublets. 

The needed primordial (B — L) asymmetry can be realized as a lepton asymmetry gener- 
ated by the out-of-equilibrium decay of heavy right-handed Majorana neutrinos, as suggested 
by Fukugita and Yanagida [Q. L is violated by Majorana masses, while the necessary CP 
violation comes about through phases in the Dirac mass matrix of the neutrinos. In a detailed 
quantitative analysis Luty showed that the scenario works for a wide range of parameters ||. 

In order to generate a lepton asymmetry of the correct order of magnitude, the right-handed 
neutrinos have to be numerous before decaying. This is only possible if they are in thermal 
equilibrium at high temperatures. In the original model |4| the right-handed neutrinos were 
only interacting through Yukawa couplings, which are far too weak to bring the neutrinos into 
equilibrium at high temperatures. Hence one had to assume an equilibrium distribution for the 
neutrinos as initial condition in previous analyses. 

An appealing way to solve this problem is to study this baryogenesis mechanism in the 
framework of an extended gauge symmetry, since right-handed neutrinos appear naturally in 
unified theories based on the gauge groups SO(10) or E 6 . As we shall see, the gauge interactions 
in which the right-handed neutrinos take part, are strong enough to bring them into thermal 
equilibrium at high temperatures. Of course, the neutrinos have to be out of equilibrium when 

x For a review and references, see m 
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decaying, i.e. the reaction rates for the gauge interactions have to fall fast enough, so that they 
cannot significantly reduce the number density of the neutrinos before they decay. 

In this paper we investigate this mechanism in the framework of an SO(10) subgroup. After 
a short discussion of the relevant Boltzmann equations in the next section, we present our 
model and calculate the needed reaction rates in section |3|. In section |] we solve the Boltzmann 
equations, first for one and then several heavy neutrino families. We explicitly show that the 
lepton asymmetry is mainly determined by the lightest of the right-handed neutrinos. Finally 
we will look at physical mass matrices for the neutrinos coming from an additional abelian 
gauged family symmetry. These mass matrices give a lepton asymmetry in the right order of 
magnitude and predict light neutrino masses and mixings of the magnitude needed to explain 
the solar neutrino deficit and a r-neutrino mass of a few eV which is needed in the cold-plus-hot 
dark matter models. 



2 Boltzmann Equations 

In a quantitative analysis of baryogenesis one can assume Maxwell-Boltzmann statistics || , so 
that the equilibrium phase space density of a particle if) with mass is given by 

(E^,T) = e~ E ^ T . (3) 

The particle density is 

n4T) = j0y 3 J d%U, (4) 

where is the number of internal degrees of freedom. The number of particles in a comoving 
volume element is given by the ratio of and the entropy density s. If the universe expands 
isentropically Y$ is not affected by the expansion of the universe, so that Y^ can only be changed 
by interactions. 

We can distinguish between elastic and inelastic scatterings. Elastic scatterings only affect 
the phase space densities of the particles but not the number densities, whilst inelastic scat- 
terings do change the number densities. If elastic scatterings do occur at a higher rate than 
inelastic scatterings we can assume kinetic equilibrium, so that the phase space density is |J 

U(E^,T) = ^e~^ T . (5) 
Consequently the Boltzmann equation describing the evolution of Y^ is (cf. || ||) 



dYj/, z x ^ 

dz sH (m^) a ^ 



JfV Y q & + a + ...^i + j + ...) 

1 ./, 'a • • • 



ye \ye q ' 7^ + j + ■ ■ ■ ~> V> + a + . • .) 



YY 2 

where z = m^/T and H (m^) is the Hubble parameter at T = m^. 
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(6) 



The right-hand side of eq. @ describes the interactions in which a ip particle takes part, 
where 7 e<? is the space time density of scatterings in thermal equilibrium. In a dilute gas we only 
have to take into account decays, two-particle scatterings and the corresponding back reactions. 
For a decay one finds 



Id ■= l eq {^ —*i + j + ...) = <?|44 f - > ( ? ) 



where Ki and K 2 are modified Bessel functions and r rs is the usual decay width in the rest 
system of the decaying particle. The ratio of the Bessel functions is a time dilatation factor. 
If we neglect CP violating effects we have the same reaction density for inverse decays, 

liD-=l eq (i+3 + — ^) = 7d • (8) 
For two body scattering one has 



y*ty + a <_> i+j + ...) = J_ / dsa^v^Ki^T) , (9) 



64tt 



where s is the squared center of mass energy and the reduced cross section a(s) for the process 
i/) + a—*i + j + ... is related to the usual total cross section er(s) by 



(jty ■ Pa) 2 ~ m^ml a(s) . (10) 



Since we have assumed kinetic equilibrium, contributions from elastic scatterings drop out 
of eq. (pf). Hence we only have to take into account inelastic processes. 



3 The model 

3.1 Gauge and Yukawa couplings 

The 16 plet of SO(10) contains, in addition to the 15 Weyl fermions of one standard model 
quark-lepton family, a right-handed neutrino Ur which is a singlet under the standard model 
gauge group. It is, therefore, natural to embed the baryogenesis mechanism of Fukugita and 
Yanagida into a SO(10) GUT. 

To explain the unification of the coupling constants in SO(10) one needs an intermediate 
breaking scale v' of the order of 10 10 to 10 13 GeV (cf. J7|). The intermediate symmetry could be 
a left-right-symmetry or a Pati-Salam-symmetry. For simplicity we take the minimal extension 
of the standard model which is based on the gauge group 

G = SU(3) C x SU(2) L x U(l) y x U(l) y , . (11) 
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This gauge group cannot explain the gauge coupling unification J?]] but it may be regarded as 
a toy model for the other symmetry groups. We then have the following breaking scheme 



SO(10) -> ... — ► SU(3) C x SU(2) L x U(l) y x U(l) y / 

{X) -^ V SU(3) C x SU(2) L x U(l)y 

^l^ v SU(3) C x U(l) em . 

Here <fi = ((p°,(p~) is the standard model Higgs doublet and x is the Higgs boson needed for 
the breaking of the U(l)y. The SU(2) L x U(l)y x U(l)y part of the lagrangian is [|| 

+ihpl L + iT7Hpv R + ieEpe R + {D^) ] {D»<j>) + (D^xV 

- (lL<f>g v VR + lL^g e eR + - x^Rhu R + h.c^j , (12) 

where we have omitted the quark fields. Il = (vl, £l) is the left-handed lepton doublet and v R 
is the right-handed neutrino. The charge conjugated field is defined by = Cv R T , where 
C is the charge conjugation matrix. g e , g v and h are the Yukawa coupling matrices which are 
responsible for the lepton masses. One can always choose the lepton fields in such a way that 
g e and h are diagonal and real. 

The covariant derivative has the form 

D^ = d,- igW, ■ f - ig'B.Y - ig's^C^Y' , (13) 

where H> 5^ and C are the SU(2) L , U(l) y and U(l)y, gauge fields. W» u , B^ and C v 
are the corresponding field strength tensors. Because of their common origin in the SO(10) 
both abelian groups have the same gauge coupling constant g' . From the structure of SO(10) 
it follows that Y' = Y — j(B — L). Therefore the fermions and the Higgs bosons carry the 
following y'-charges 

y' (h) = \ V (e R ) = \ Y' (u R ) = \ 

Y'(q L ) = -\ Y'{u R ) = \ Y'(d R ) = -\ (14) 

Y'{*) = -\ Y'( X ) = - 5 -. 
The most general Higgs potential is 

V( X , 0) = /ii 0V + //a X ] X + \\x (0V) 2 + (x ] xf + A 3 (x ] x) (j> ] <P) , (15) 

where A 3 > — \f\\\2, so that the potential is bounded from below. Spontaneous symmetry 
breaking leads to two massive neutral gauge bosons Z and Z' . 



The right-handed neutrinos acquire Majorana masses M = hv' when the U(l)y is broken. 
Additionally the neutrinos get Dirac masses mo = vg u with v = 174 GeV at the electroweak 
symmetry breaking. This offers the possibility to explain the smallness of the u e , and v T 
masses via the seesaw mechanism [[J. 

Since we want to explain the generation of a lepton asymmetry before the electroweak 
phase transition, we can take ((f)) — v — 0, so that all mixing angles vanish. Therefore, all the 
standard model particles are massless and the additional neutral gauge boson Z' is identical to 
the vector field C M . For (x) — v' ^ the mass of the Z' is 

5 

Since Yukawa couplings are usually small when compared to gauge couplings, one expects that 
the Majorana masses Mj of the right-handed neutrinos are small compared to the Z' mass. 
The Higgs boson H f contained in the field x gets the mass 



m H ' = yX 2 v' . (17) 



To realize the assumed hierarchy v' ^> v one needs A2 > 2 ■ 10 3 ||10|| . Hence m#£ is of the order 



of mz'i so that processes in which a H' Q takes part are kinematically suppressed at temperatures 
below uiz'- Moreover the Yukawa couplings h of H' Q were assumed to be small in comparison 
with the gauge couplings. Therefore the Higgs boson H' Q can be neglected in the following. 

Since we have assumed that the matrix h is real and diagonal the weak eigenstates are equal 
to the Majorana mass eigenstates, 

v = v L + v C L and TV = v R + v% , (18) 

where v and iV are four-component Majorana spinors. 

3.2 Reaction rates for lepton number violating processes 

We are now able to calculate the relevant reaction rates. It appears reasonable to assume 
a mass hierarchy of the form M\ <C M2 -C M3 for the right-handed neutrinos, so that the 
lightest right-handed neutrino, iV 1 , will still be in equilibrium when iV 2 and iV 3 decay. The 
lepton number violating interactions mediated by iV 1 will, therefore, wash out an asymmetry 
generated by the decays of the neutrinos iV 2 and iV 3 . Hence a significant lepton asymmetry can 
only be generated by A^ 1 -decays and we expect that the temperature at which the asymmetry 
is generated is of the order of M\. Consequently it is useful to relate all the masses and energies 
to M\. We define the following dimensionless quantities 



M ^ v-^' and t " — S riQl 

The lepton asymmetry is generated by the CP violating decay of the right-handed neutrinos. 
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Figure 1: Right-handed neutrino decay at tree and one-loop level. 



Making use of the relation 



1 9 

Qv = m D - = m D 

v V2M W 



(20) 



where Mw = 80 GeV is the W boson mass and g is the SU(2) coupling constant, one finds for 
the decay width at tree level H 



:= f „ (n 3 



+ l)+ T rs (N j 



+ 1 



a Mj (mJ D m£)) 



33 



sin 2 9 4 



M 2 W 



(21) 



Here 9 is the weak mixing angle. The leading contribution to the CP-asymmetry in the decay 
of N 3 reads 

f „ (n 3 -> 0t + /) _ f rs (nj + T 



e 5 :-- 



(22) 



f„ (Ni -> 0t + z) + f„ (iVi -^0 + 1) ' 

It is due to the interference between the tree level decay amplitude and the one loop amplitude 
shown in fig. [TJ. One finds^] 

- — — y^Im (m^mo) 
'1 + x" 



a 



4M 2 v sin 2 9( m ^m D ) :jj ^ 



(23) 



with f(x) = \fx 



1 - (1 + x) In 



(24) 



From the definition ( |2"2"D of the CP asymmetry it follows that the reaction densities for decays 
and inverse decays can be parametrized in the following way 



with 



7 e<? (^3 _^ 0t + (j = + J _> ^ 

7 e<? (iW -> + 7) = 7 eg (0 f + / -> N 3 ) 

lDj : = rfl (N3 _ 0t + ^ + _ ^ + fj 



-(l-e J ) 7Dj 



Dj ^ K 2 (M,/T) 



(25) 
(26) 
(27) 



Since the CP-asymmetry is a higher order effect we have to take into account other higher 
order effects as well, especially the L violating lepton Higgs scattering I + — > I + <ft shown in 



The results quoted in the literature differ from our r esult by a factor 8 (ref. §), 2 (ref. 0) or 18 (ref. {§) 
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Figure 2: Lepton number violating lepton Higgs scattering 



fig. |2| with a right-handed neutrino as intermediate state. As shown in ref. |J this diagram is 
necessary to avoid the generation of a lepton asymmetry in thermal equilibrium. 

However, the scatterings with a real intermediate neutrino were already taken into account, 
because they can be described by an inverse decay followed by a decay. Therefore, we have to 
remove the contributions from physical intermediate states from the reduced cross section|J. 

In H the flavour structure of this diagram was neglected, i.e. only the contribution from 
the lightest right-handed neutrino iV 1 was taken into account. We will take into consideration 
two right-handed neutrinos, N 1 and N 2 . When discussing our results we will see that the third 
neutrino can indeed be neglected. The reduced cross section is then given by: 



2tt 1 



a 



~ \J2 a 3 ( m o m o J 



\ 2 
33 



X 

— + 



2x 



+ 



x 



2D){x) 



1 + 2 



x + a 



3 



In 



x + a. 



+2y / aia2 Re 
(x + ai)(x + ai 



(rn} D m D 



x 



+ 



X 



+ 



X 



DAx) D 2 (x) 2D 1 (x)D 2 (x 



2a 2 ) . f x + a\ 
- In 



(x + a 2 )(x + a 2 - 2ai) fx + a 2 

- In 



• (28) 



D 2 (x)(a l -a 2 ) "\ at J D 1 (x)(a 2 -ai) \ a 2 

The prime denotes that we have subtracted the contributions from real intermediate states and 



1 



Djix) 



(x — Oj) 2 + QjCj 



with Cj :-- 



1 Dj 

Mi 



(29) 



is the off-shell part of the propagator. 

There are some other L violating processes mediated by a right-handed neutrino in the 
t-channel, like the scattering I + I — > <p + <fi. The reduced cross section for this process, which 
was neglected in M, is 



o-nAs) 



+ Re 



\ 2 

33 



X 



1 



x + aj x + 2aj 



In 



X + OLa 



a-. 



(30) 



(™>D 



m D 
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y/aia 2 



(ai - a 2 )(x + ai + a 2 ) 
The same result is valid for the process I + 1 



' x + 2ai)\n[^^] -(x + 2a 2 )ln / '" ' 



+ 0t and the back reactions. 
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t > 1 > Q 



Figure 3: Lepton number violating scatterings mediated by a standard model Higgs boson in the 
s- or t-channel 

The Yukawa couplings of the standard model are very small and can be neglected, with 
the exception of the term responsible for the top quark mass m t . Thus we have to take into 
account the L violating interactions of a right-handed neutrino with a top-quark. There are 
s- and t-channel contributions, shown in fig. || The reduced cross section for the s-channel 
process N J ' + I — > t + q is || 



One finds the same result for the process iV- 7 + 7 — > t + q and the corresponding back reactions. 
To get a well defined result for the t-channel contribution one has to introduce a Higgs mass 
Keeping only the leading contributions for —>■ one finds [|| 



„ i , . 37ra 2 m? 

oi + (s^ - 



M D m D ) . . 



H — - In ' 



(32) 



x x \ y' 

with the dimensionless squared Higgs mass y' := m^/Mf. Since these processes have only a 
small influence on the generated asymmetry, the results are very insensitive to the value of m^. 
In the following we use m<$, = 800 GeV. 

3.3 Reaction rates for lepton number conserving processes 

The lepton number violating processes are too weak to bring the right-handed neutrinos into 
equilibrium at high temperatures. Therefore, one has to consider lepton number conserving 
processes, which were neglected up to now, like the iV- 7 pair creation and pair annihilation. The 
dominating contribution comes from the Z' exchange in the s-channel shown in fig. |]a. The 
total reduced cross section for the processes / + / <-> + , where / is a standard model 
fermion, and <fi + (jr <-> + is 

A 4225vr a 2 y/x 3/2 

CZ'(S) — 7-— -. rr: [X — AcLs) ' . 33 

K ' 216 cos 4 # (x-y) 2 + yc v Jj K ' 

To get a well defined result we have introduced the Z' decay width into the propagator, 

(34) 




~ am z , 
with T z > = — 57: 

COS^ V 





4a, V „, ,,169 
- 4 ^ )+ 144 




Figure 4: Lepton number conserving processes: 

(a) pair creation 

(b) process N* + N* -> iV J ' + JV 3 ' (i ^ j) 
(c,d) eZasizc scattering N { + N* -> N* + N { 



Since we will have different heavy neutrino flavours we have to consider transitions between 
right-handed neutrinos like iV^ + iV* <-» + (i ^ j) with an intermediate Z' in the s-channel, 
shown in fig. £|b. The reduced cross section for this process is 



„ , 6257T a 2 1 (x — 4a,) 3 / 2 (x — 4a,) 3 / 2 a,-a,- /x — 4a,- /x — 4a,- . 
216 cos 4 ft i (sc — yj^ + yc y z y x y x 

Finally we have to look at the elastic processes to check that the assumed kinetic equilibrium 
is a good approximation. For this purpose we have considered the process N l + N l — > N l + N l . 
The dominating contribution to this transition comes from Z' exchange in the s-channel (fig. |]c) 
and in the t- and w-channel (fig. f|d). The reduced cross section for this process reads 

, , 625vr a 2 [ 1 {x-Aaif 

O-el(s) 



72 cos 4 6 { 3x (x — y) 2 + yc 

| x - 4a^ y(x - 4a^) 2 + 2(y - 2^) 3 + (8a 2 + 3y 2 )x - Aa^y 2 + yoj + 4a 2 ) | 
sy 2 x — 4aj + y 

| 1 (3y-4a l )(x-4a t ) 2 + (3x-20a l )y 2 + 2y(y 2 + 8a 2 ) ^f y \\ 
xy x — Aai + 2y \x — 4a, + y J J 



4 Results 
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4.1 Constraints on the parameters 

Before solving the Boltzmann equations we will try to constrain the parameters by looking at 
the reaction rates. The neutrinos have to be out of equilibrium when decaying, i.e. they have 
to decouple before decaying. Therefore, the decay rate To has to be smaller than the Hubble 



parameter H at temperatures T « Mj. This gives the following constraint [11, 12 



T D (T = Mi) < 3H(T = Mi) <£> m x := ( m ^)ii < g _ 1Q _ 3 eV _ (3?) 

rh\ is the mass of the lightest neutrino, if mjj is approximately diagonal. 

Furthermore the pair annihilation rate IV for the neutrinos has to be smaller than 3H at 
T w Mi. This gives a lower bound on the Z' mass 0, [H| 



/ M \ 3/4 

r z ,(T = Mi) < 3H(T = Mi) ^ m z , > 2 • 10 11 GeV f ^L— J • (38) 

On the other hand the L asymmetry must be generated before the electroweak phase transition, 
because we need the sphalerons which are no longer in equilibrium at temperatures below 
100 GeV. Since the electroweak phase transition takes place 10~ 12 s after the big bang, the 
lifetime of the right-handed neutrinos has to be smaller than 10~ 12 s. This yields the following 
constraint 

(1 n 10 PpV\ 
■ m 

Finally the reaction rates and T^t for the lepton number violating scatterings mediated 
by the right-handed neutrinos must not wash out the generated lepton asymmetry at low 



temperatures. This gives the following condition |TB 



(1 10 CYV\ 
w ' ) • (40) 

The sum over m 2 can be interpreted as sum over the squared masses of the light neutrinos, if 
mo is approximately diagonal. 

We have shown the different reaction rates and the Hubble parameter in fig. |5|, where we 
have chosen the following parameters 

fh\ = 10~ 3 eV , Mi = 10 10 GeV and m z > = 2- 10 11 GeV . (41) 

The reaction rates have the correct behaviour. Tz 1 is much bigger than 3H at temperatures 
T ^> Mi, so that the neutrinos iV 1 come into equilibrium at high temperatures, while Tz> 
is small enough at T w Mi to avoid the lepton number conserving pair annihilation of the 
neutrinos. Furthermore one can see that the reaction rates taken into account in previous 
analyses [|], i.e. T^, Tjd, Tjv, T^ )S and are far too low to bring the neutrinos into 
thermal equilibrium at high temperatures. 
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z=M 1 /T z = M,/T 

Figure 5: Comparison of different reaction rates with the Hubble parameter H. T z> is the pair 
annihilation rate, the decay rate, Tju the reaction rate for inverse decays, and T^ t for 
lepton number violating processes with an intermediate N 1 , and T^t for the interaction of 
an incoming N 1 with a top-quark and T e i for the elastic N 1 -scattering. 

Since the decay rate is smaller than 3H for T w Mi, the neutrinos have enough time 
to deviate from the thermal equilibrium, so that the out-of-equilibrium condition is fulfilled. 
Moreover the reaction rates Tjy and for the lepton number violating interactions and IYd 
for inverse decays are so small that the L asymmetry cannot be erased at temperatures T < Mi. 

Finally we have plotted the reaction rate T e i for the elastic scatterings N 1 + N 1 <-> iV 1 + N 1 . 
As one can see T e i is much bigger than Y z> in the relevant temperature region T rs Mi, so that 
one can assume kinetic equilibrium. 



4.2 Decay of the lightest heavy neutrino 

For simplicity we begin by neglecting the decays of the right-handed neutrinos iV 2 and iV 3 . 
Using the formalism of section ^| and the reaction rates calculated in section ^| we can immedi- 
ately write down the Boltzmann equations. If we assume that all the standard model particles 
are in thermal equilibrium, the Boltzmann equations for the iV 1 number and for the (B — L) 
asymmetry are 





- * J 




dz 


sH{Mi) 




dYB-L 








d^ 


sH (Mi) I 




2 Y t eq 




Jn 1 



(42) 
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Figure 6: Typical solutions of the Boltzmann equations. 



Yi 



B-L 



2 lN + 2 lN , t + 2 7 J_ t + 



N 1 



(43) 



These equations have to be solved numerically. We start at T ^ Mi with the initial conditions 



Y N i — Yffj, 



and 



B—L 







(44) 



and follow the evolution of these quantities down to low temperatures. We want to emphasize 
that the results are independent of the initial conditions (f44f) , because the reaction rates are 
so high that the neutrinos are rapidly driven into equilibrium at high temperatures while any 
primordial (B — L) asymmetry is washed out. Typical results are shown in fig. || where we 
have used the following parameters 

Mi = 10 10 GeV m z > = 10 13 GeV 



and 



nil 



(m^m D ) 



ii 



Mi 




ei = -5-10' 

for fig. 6a 
for fig. ^3 . 



(45) 
(46) 



One recognizes very clearly the departure from thermal equilibrium and the generation of the 
asymmetry in the correct order of magnitude. If the neutrinos decay at very low temperatures 
T ^ Mi, the back reactions and the lepton number violating scatterings can be neglected, so 
that we expect a (B — L) asymmetry of the order |J 

Y B -l « SiY m (T = Mi) = 2.1 ■ 1(T 10 
This is exactly the result of fig. |]a, while one has 

Yb—l = 1.6 ■ 10" 
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(47) 



(48) 
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Figure 7: The generated (B — L) asymmetry for mz> = 10 3 Mi (a) and m z > = lOMi (7>J and for 
Mi = 10 s GeF (dotted line), M 1 = IO 10 GeV (solid line) and Mi = 10 12 GeV (dashed line). 

in fig. |6|b. This means that the inverse decays and the lepton number violating scatterings 
already erase a part of the asymmetry. A further increase in {Trv D rriD)ii would intensify the 
dissipative processes, so that the final asymmetry would be even lower. 

One gets almost the same results if one varies the parameters and keeps the ratios y = 
(mz> /Mi) 2 and fh\ and the CP asymmetry E\ constant. In fig. |^a we have plotted the generated 
asymmetry as a function of rh\ for y = 10 6 and E\ = —5 • 10~ 8 for the neutrino masses 
M x = 10 12 GeV, 10 10 GeV and 10 8 GeV. 

First of all one notices that the condition Q57D delimitates quite accurately the range of rh\ 
in which the generation of the (B — L) asymmetry is possible with a reasonable choice for the 
CP asymmetry. 

In the second place it is remarkable that the asymmetry depends almost entirely on rh\ and 
not on {jrv D mj))\\ and M\ separately. This becomes clear if one looks at the right-hand side of 
the Boltzmann equation fl43l) . Indeed the decay term depends only on rhi, 

z 



7W(M7) lmCCmi - 

The same proportionality holds for 7^.^ and 7^. On the other hand one finds 



• 7at oc m 2 Mi 



and 



■ Tjv t oc m\ Mi 



(49) 



(50) 



si? (MO <" - 1 ' ^ - sH{M l ) 

for the lepton number violating scatterings. Since these terms can erase a lepton asymmetry 
one would expect that the asymmetry is falling with M\. This is exactly what one can observe 
in fig. [7|a. 
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Up to now we have always chosen a large Z' mass, mz> = lO^Mx, to ensure that pair 
annihilation processes do not influence the evolution of Y^i at temperatures T < M\. As a 
counter-example we have repeated our calculations with a lower Z' mass, mz> = lOMx- The 
results are summarized in fig. 0b. In contrast to fig. 0a one recognizes a strong dependence of the 
asymmetry on M\. This comes about because of the pair annihilation term in the Boltzmann 
equation (|42), which is inversely proportional to M±, 

z 1 



sH (MO" 7 * K Ml ■ (51) 

Therefore, the pair annihilation processes should be more effective at low neutrino masses M\. 
This explains the dependence of the curves in fig. 0b on M\. 

Fig. 0b may further be used to confirm the lower bound ( p8|) on the Z' mass. One can 
generate the required asymmetry for M\ = 10 12 and 10 10 GeV but not for Mi = 10 8 GeV, 
because the condition (|38| ) is only fulfilled for the first two values of Ml. 

The curve for Mi = 10 8 GeV in fig. 0b has another interesting property. Up to now the 
generated asymmetry was always monotonically falling with fhi, while Yb-l has now a max- 
imum for fhi ~ 10~ 3 eV. The reason is, that for this value of fhi the neutrinos decay before 
they can be depleted by pair annihilation processes, while the reaction rate for inverse decays 
is still too small to erase the asymmetry. The evolution of the asymmetry is hampered by pair 
annihilation processes at lower values of fhi and by dissipating processes like inverse decays at 
higher values of rh\. 



4.3 Decays of two heavy neutrinos 



Now we can refine our results by adding a second family of right-handed neutrinos. The 
Boltzmann equations are straightforward generalizations of the equations ( [42]) and (|43|) 

2 



dY } 



B-L 



dz 



1 Yb-l Y Ni 
+ £ j ' 



2 Y t eq 



1 No 



i Yb-l ro , o i i Yb-l 



1 i 



E 



\req 



Y N j 



1 Ni 



(52) 



dY N i 
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4>,t 
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Y, 



N 1 



Jn 1 



lz> + 



Y, 



N 1 



Jn 1 



'Y^ 



7AT1AT2 



(53) 



dY N 2 
dz 
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— 1 



7 D2 + 2 7 J, S + 4 7 J i4 
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* TV 2 , 



72' + 
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» TV 2 , 



Yi 



N 1 



7AT1JV 2 



(54) 



(a) 



(b) 




where Jn 1 ^ is the reaction density for the ./V^-iV 2 -transitions described by the reduced cross 
section (|35|) . 

The constraints which were derived in section |4.1| remain valid, but the parameters for the 
two heavy neutrinos are not always independent of each other. For the CP asymmetries one 
has for example: 



Ex oc Im 



{m l D m D 



and e 2 oc Im 



(55) 



Since 



M^M, , (56) 

the two CP asymmetries have a different sign. Therefore, the asymmetries generated in the 
decays of N l and N 2 will also have a different sign, so that they can cancel each other if 
the neutrinos are mass degenerate. But this should not be a problem if the neutrinos have 
hierarchical masses like the charged fermions, because the asymmetry generated by the heavier 
neutrinos is more affected by the dissipative processes and because the mass dependence of e\ 
and E2 ensures that \e^\ is smaller than |ei| if M2 > M\. Analogous phenomena have already 
been observed in other models fil l. 

Typical results of the numerical integration of the Boltzmann equations are shown in fig. || 
In both cases we have chosen 



Mi = 10 10 GeV 



Mo = 10 11 GeV 



M z > = 10 13 GeV and Re 



\ 2 



m J D m D j 



12 



0. 

(57) 
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The results are almost independent of the choice Re (m^rriD^ = 0, because this parameter 
only appears in interference terms between the different diagrams contributing to the lepton 
Higgs scatterings which have only a small influence on the generated asymmetry. 
Furthermore we have for fig. [8k 



m 2 



Im 



m x = — = l(T b eV 



100 
£l = -5 • 10" 



and 
and 



J 2 



Ml 
e 2 = 3.7-10 



7.7 ■ 10~ 8 (eV) s 



-10 



(58) 



and for fig. Rib we have 



mi 



m 2 

Too 



Im 



10~ 4 eV and 



(m} D rriD) 



12 



& e x = -5 • 10 



-8 



and 



Mi 

e 2 = 3.7- 10 



7.7 • 10~ 6 (eV) s 



-10 



(59) 



Therefore, we have chosen exactly the same parameters for N l as in the one-family calculations 
of fig. [], so that we can compare the results. First of all one sees that the asymmetry changes 
sign at z ~ 6 and z ~ 0.8 as we had predicted it at the beginning of this section. But the 
influence of the second neutrino family on the final asymmetry is only very small and we recover 
the results (|47|) and (fE|). Consequently one can always neglect the heavier neutrino families if 
one has a pronounced mass hierarchy. 



4.4 Physical mass matrices 

Up to now we have only verified that the parameters of the theory can be chosen in such a way, 
that they allow the generation of the requested asymmetry. However, in a physical theory the 
parameters have to comply with some other conditions. First of all the mass scale of the right- 
handed neutrinos and the Z' boson have to be related to the intermediate breaking scale v' . 
Next the predicted light neutrino masses have to be consistent with the experimental bounds. 
We will focus on the following conditions: 

1. The solar neutrino deficit can be explained by the Mikheyev-Smirnov-Wolfenstein (MSW) 
effect, if the v e and masses and their mixing angle fulfill the following conditions (cf. 0) 

Sm 2 := m 2 ^ - m\ = (0.3 - 1) • 10~ 5 (eV) 2 , (60) 
sin 2 2# eM = 0.003 - 0.012 ; (61) 

2. a v T mass of a few eV is needed in the cold-plus- hot dark matter models. 

A realistic pattern of masses and mixings of the known fermions can be obtained based on a 
gauged abelian family symmetry broken below the unification scale This model predicts 
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the following mass matrices for the up and down quarks (suppressing unknown phases and 
factors which are assumed to be of order 1) 



M„ 



/ e |2+6a| £ 

e |3a| 
I £ |l+3a| 



3« 
.2 



|l+3a| 

e 
1 



/ g|2+6a| ^|3a] g|l+3tt| \ 



m t 



and 



J 



p|3a 
|l+3o| 



(62) 



/ 



with a = as/ («2 — ^3); where a, is the family charge of the z-th generation quarks. One obtains 
a good agreement with the measured values of the quark masses and mixings for 



1 



and 



0.23 . 



(63) 



Similarly the leptons of the z-th generation have the family charge eij, and therefore the mass 
matrix of the charged leptons reads 



/ gl8-»| ^3 



Mi 



e 3 

|4-6j 



e|4-6| \ 

g|2(l-b)| 



for integer b 



1 J 



«2 — «2 

a 2 — a 3 



(64) 



The known masses of the charged leptons can be explained for 6 = 0. Mi is diagonalized by an 
orthogonal matrix Ri. If one neglects higher order terms in e, Ri is given by 



Ri. 



6, 



0/Lf 



e/j. 
1 



O (e 4 ) \ 
O(e) 

1 / 



with 5 eM = J — = 0{?) . 



(65) 



This model also predicts Dirac masses mo for the neutrinos and Majorana masses M for the 
right-handed neutrinos [[Ljj, 



m D 



( e 16 e 6 
e 6 e 4 
\ e 8 ere 2 
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( e 8 re 3 e 4 \ 
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P I 



f3 Af2 
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re 3 Ce" ' 



M 3 , 



(66) 



where we have introduced four factors p, a, r and ( of order 1, which will be used to fix the 
physical parameters. As we are working with mass eigenstates N l we have to transform mjj 
and M into a basis in which M is diagonal. Up to higher order terms in e the matrix Rm which 
diagonalizes M is 



/ 1 



R 



M 



e 2 e 4 \ 



— e 



1 e 

which yields the following masses for the heavy neutrinos 
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Figure 9: Solutions of the Boltzmann equations with physical mass matrices. The equilibrium 
distributions are represented by dashed lines. 



Therefore the mass ratios of the right-handed neutrinos are 



(M 2 



(C 



4r 2 e« 

The Dirac mass matrix in the new basis is 

( 



and 



«3 



M 
Mi 



4r 2 e 10 



m' D := rriD R 



M 
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e 6 
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V -(p + a)e 4 -(p-cre)e p ) 



The mass parameters which appear in the decay widths are thus given by 
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m'^ D m' D 
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(p + afe 2 m 2 u:j 



Mi 
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M 3 
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(p - ae) 2 ml 3 



M, 



|C-1| M 3 



(69) 



(70) 



(71) 



Up to now we have neglected all the phase factors that we need for CP violation. In a Majorana 
mass matrix the number of physical phases for the n generation case is n(n — l)/2 (cf. 0). 
They can be chosen in such a way, that the squared non diagonal elements of m D m' D are purely 
imaginary. Then the CP asymmetries are 



Si ps r 



IC 



^2 



"3 



3.6 • 10 4 GeV 



Now we have to specify the parameters p, a, r and £. A good choice is 



p = 0.9 



a 



t = —— and C = 2 . 
4 



(73) 
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The mass scales are approximately fixed by their breaking scale, 

m U3 = 50 GeV , M 3 = 5 • 10 11 GeV and M z , = 10 12 GeV 
Then the Majorana neutrinos iV 1 and N 2 have the masses 

Mi = 3.1 • 10 8 GeV and M 2 = 2.6 • 10 10 GeV . 
For the coupling parameters ( [FID and the CP asymmetries (|72[) one finds 

mi = 8 • lCT 4 eV and m 2 = 6.4 eV , 
ei = -1.3 • 10~ 6 and e 2 = 4.6 • 10~ 5 . 



(74) 

(75) 

(76) 
(77) 



The corresponding solutions of the Boltzmann equations are shown in fig. |^. The generated 
asymmetry, 

Y B _ L = 5 ■ 1(T 10 , (78) 

is even bigger than requested, but it is always possible to reduce Yb-l by choosing the phase 
factors in an appropriate way. One finds similar results for other values of the parameters p, 
cr, r and (. 

Finally we have to check if these matrices can predict the desired light neutrino masses and 
mixings. The Majorana masses of the light neutrinos are given by the seesaw formula [0], 
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-2e 15 (a - 1)^ 
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2 —4 

a t e par e 



pe 7 \ 
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(79) 



The eigenvalues of this matrix are 
8.9- 10~ 7 eV , 



m,. 



1.9-10 _3 eV and m h 



2.2 eV 



(80) 



2e 3 
1 



e 5 



These values fulfill the MSW condition (|60| ) and the v T mass is in the correct range for the 
cold-plus-hot dark matter models. m u is approximately diagonalized by the following matrix, 

/ 1 

Ry w -2e 3 

V -6 5 -(a/p)f 1 

The leptonic analogon of the Cabbibo-Kobayashi-Maskawa (CKM) matrix is 



-2 



(81) 
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5 etl - 2e 3 
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(82) 



Therefore the mixing angle complies with the MSW condition 



sin 2 29, 



e/j, 



0.009 



(83) 
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5 Conclusions 



We have seen that the cosmic baryon asymmetry can be explained by the lepton number 
violating decays of heavy Majorana neutrinos combined with the anomalous electroweak (B+L) 
violation. The lepton asymmetry generated in our model is independent of the initial conditions 
on the heavy neutrino density, which was not the case in a previous analysis of this mechanism 
. This is a consequence of the new gauge interaction which we have introduced, and which 
is related to the spontaneous breaking of lepton number. 

We have also considered the decay of more than one heavy neutrino and have seen that 
the generated asymmetry is determined by the properties of the lightest heavy neutrino, if the 
right-handed neutrinos have a pronounced mass hierarchy. 

By performing explicit calculations we could show that the generation of a lepton asymmetry 
is possible at every temperature between the intermediate breaking scale and the electroweak 
breaking scale. Furthermore we have checked that neutrino mass matrices may explain the 
generation of the asymmetry and low energy neutrino phenomenology at the same time. 

In supersymmetric models of inflation the reheating temperature has to be lower than 10 5 to 
10 8 GeV to solve the gravitino problem [17], [T3|. Therefore in these models baryogenesis has to 
take place at relatively low temperatures. Since this is possible in our model, a supersymmetric 
generalization should be viable. 
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